On the Ejection Mechanism of Bullets in SS 433 
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ABSTRACT 



We discuss plausible mechanisms to produce bullet-like ejecta from the pre- 
cessing disk in the SS 433 system. We show that non-steady shocks in the sub- 
Keplerian accretion flow can provide the basic timescale of the ejection interval 
while the magnetic rubber-band effect of the toroidal flux tubes in this disk can 
yield flaring events. 

Subject headings: accretion, accretion disks — hydrodynamics — instabilities — 
shock waves — stars: individual (SS 433) — stars: mass loss 
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1. Introduction 



SS 433 remains one of the most enigmatic objects in the sky. Even twenty-five years 
after its first appearance in the catalogue of Stephanson & Sanduleak (1977), it is not clear 
whether the compact object is a black hole or a neutron star. However, there is ample 
evidence that the companion is an OB type star with an orbital period of 13.1 days, which 
is losing mass at the rate of about lO'^M© yr~^ (van den Heuvel 1981), corresponding to 
extremely super-Eddington accretion regardless of the mass of the compact object. 

One of the most curious properties of the jets of SS 433, which first made their presence 
distinctly felt through the emission of variable Ha lines, is that they are apparently ejected 
as bullets (e.g. Borisov & Fabrika 1987; Vermeulen et al. 1993; Paragi et al. 1999; 2002; 
Gies et al. 2002), with a surprisingly nearly constant radial velocity of about 0.26c. The 



- 2 - 



absence of a significant intrinsic rotational velocity (i.e., v^j)) component is clear from the 
fact that the kinematic model (e.g. Abell and Margon 1979), which assumes only radial 
injection, quite accurately explains the time variation of the red- and blue-shifts of the Ha 
emission from the jets with a period of 162 days, which is attributed to the precession of 
the accretion disk about the compact object. The radial velocity is less than the maximum 
allowed sound speed of c/ y/S and thus hydrodynamic acceleration could, in principle, explain 
it. Therefore one may not require a magnetic or electrodynamic acceleration processes (e.g. 
Belcher & MacGregor 1976; Lovelace 1976). However, the rather good coUimation (Margon 
1984; Paragi et al. 1999) supports the hypothesis that a substantial degree of confinement 
produced by toroidal flux tubes may be present. Gies et al. (2002) showed that the ratios 
of the Ha emission equivalent widths from the approaching and receding jets as a function 
of precessional phase only could be nicely fit if these emission components are bullet-like. 
Indeed, the recent Chandra X-ray Observatory discovery of X-rays at a distance of about 
lO^^cm from the center may result from the collision of such bullets (Migliari, Fender, & 
Mendez 2002). 

SS 433 poses another interesting problem: it was pointed out by Chakrabarti (1999) and 
Das & Chakrabarti (1999) that significant outfiows are produced only when the accretion rate 
is such that the X-ray source is in a low/hard state, and all the observational indications 
in other micro-quasars also suggest that the jets are indeed produced in low/hard states 
(Corbel et al. 2001; Klein- Wolt et al. 2002). However, it is difficult to imagine how SS 433 
manages to remain in the low/hard state with 10~^Mq yr~^ of wind matter ejected from 
its companion. The answer to this quandary probably lies in the recent results of Paragi et 
al. (1999) and Blundell et al. (2000), whose high resolution radio maps show that there is a 
large region of roughly 50AU in radius which is filled with enough gas and dust to obscure 
the accretion disk and the base of the jets. They also found an equatorial outfiow. Gies et 
al. (2002) present additional evidence from observations of the "stationary" Ha and He I 
lines for an extended "disk wind'' . So it is distinctly possible that most of the matter from 
the donor is rejected either by centrifugal force (Chakrabarti 2002) or by radiation force far 
outside the central accretion disk, and thus the compact object receives only a few times the 
Eddington rate (Meci) of its companion's wind matter to accrete. This consideration finds 
further support from the fact that the kinematic luminosity of the jet itself is around 10^^ 
erg s^^ (Margon 1984), which corresponds to about one Eddington rate for a IOMq compact 
object. 

In numerical simulations of supercritical winds by Eggum, Coroniti & Katz (1985) 
designed to model SS 433, it was shown that only a fraction of a percent of the infalling 
matter is ejected from a radiation pressure supported Keplerian disk, which indicates that the 
accretion rate must be at least IOOMeci if the accretion takes place through a Keplerian disk. 
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On the other hand, numerical simulations of a sub-Keplerian disk by Molteni, Lanzafame & 
Chakrabarti (1994) suggest that about 15 — 20 percent of matter is ejected as an outflow, 
indicating that the accretion rate onto the compact object in SS 433 need be at most a 
few Mgd- Similar simulations with different parameters yield situations where no steady 
shocks can form, even though two saddle-type sonic points arc present (Ryu, Chakrabarti, & 
Molteni 1997, hereafter RCM); under these conditions large scale shock oscillations produce 
intermittent outflows instead of continuous outflows. Since the compact object is a wind 
accretor, a low-angular momentum, sub-Keplerian flow is the most likely description of the 
accretion flow. Indeed, the presence of sub-Keplerian flows in several other high mass X-ray 
binaries has now been verifled (Smith, Heindl & Swank 2002). 

In this Letter, we present a few scenarios leading to ejection of matter as bullets in SS 
433. We discuss four possible ways to create blobs of matter emerging from the disk and 
conclude that periodic ejection of the blobs by the large scale oscillation of an accretion shock 
(something like a piston) may be the fundamental production mechanism of the "normal" 
bullets. The irregularly observed rapid flaring (Vermeulcn ct al. 1993) could be understood 
in terms of the catastrophic collapse of toroidal magnetic flux tubes, very similar to what has 
been argued to be occurring in GRS 1915+105 (Vadawale et al. 2001; Nandi et al. 2001). In 
the next section we discuss these processes and their suitability or unsuitability for SS 433. 
In §3, we present concluding remarks. 



2. Mechanisms to produce buUet-Iike ejecta from accretion flows 

In both the works of Eggum et al. (1985) and Molteni et al. (1994) continuous ejection 
was reported when a radiation pressure dominated Keplerian disk, or a sub-Kcplcrian disk 
capable of producing a steady shock, were considered. However, in SS 433 the basic ejection 
is buUet-hke and since the size of the X-ray emitting region is smaller than ~ 10^^ cm 
within which the material in the jets is already accelerated to Vjet ~ 0.26c (Watson et 
al. 1986; Stewart et al. 1987), the bullets are not expected to be delayed by more than 
Ix/vjet ~ 100 s. Indeed, recent Rossi X-Ray Timing Experiment (RXTE) observations of 
hard X-rays from SS 433 indicated variability on time scales of 50-1000 s (Safl-Harb & 
Kotani 2002), roughly corroborating this picture. In fact, a simultaneous measurement of 
a flare at 2 GHz in the radio (Kotani & Trushkin 2001) and in hard X-rays (Safl-Harb & 
Kotani 2002) indicated a strong anti-correlation of radio and X-ray fluxes, similar to what 
is observed in GRS 1915+105 (Mirabel & Rodriguez 1994). Moreover, the X-ray luminosity 
is very low (~ lO^^erg s~^) and is believed to come from the base of the jets (Watson et al. 
1986). It is beheved to have a thermal origin and EXOSAT (Watson et al. 1986) and GINGA 
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(Yuan et al. 1995) observations were adequately fitted with a thermal bremsstrahlung model 
with kT ^ 30 keV. The overall spectral shape suggests that the source has always been in 
a standard low/hard state and so far no quasi-thermal emission expected from a 'Keplerian 
disk' has been detected. From the interaction of the jet with the supernovae remnant W50, 
the lower limit of kinematic luminosity is found to be at least lO^^erg (Biretta et al. 
1983; Davidson & McCray 1980). This means that the mean mass outflow rate is around 
lO^^g s~\ and if most of it is in the form of bullets ejected at 50-1000 s intervals, the mass 
accumulated in each bullet should be in the range of 10^^-10^^ g. 

The above data implies that the essential features that one must explain when attempt- 
ing to produce bullets out of the accretion disks are: (a) the disk should be a sub- Keplerian 
flow; (b) the object (black hole or a neutron star) and its surroundings should be in a 
low/hard state; (c) bullets should be ejected in 50-1000 s time-scales under normal circum- 
stances; (d) the mass of each bullet should be around lO^^^^^g; and finally, (e) there should 
be occasional flaring with an anti-correlation of radio and X-ray emission. We now discuss 
several scenarios and present what we believe to be the most probable picture of what is 
going on in SS 433. The four processes are schematically shown in Fig. l(a-d). 

2.1. Cooling of the jet-base by Comptonization and separation of blobs 

It was shown by several numerical simulations that significant outflows are produced 
from regions very close to the inner edge of the accretion flow, possibly from the centrifugal 
pressure dominated region (Molteni et al. 1994; Molteni et al. 2001). These jets are launched 
subsonically, but quickly pass through the inner sonic point to become supersonic. In the 
subsonic region while the matter moves slowly, the density is high and the optical depth 
could be large enough (r > 1) to undergo Compton cooling (Fig. la) provided there is a 
Keplerian disk underneath to supply soft photons. A part of the outflow, which was subsonic 
previously becomes supersonic because of this rapid cooling and separates from the base 
of the jet. This separation of blobs is expected to occur at the sonic surface Tc which is 
~ 2 — 3rs, where, is the size of the centrifugal barrier; see Chakrabarti (1999). 

This possibility, though attractive, and in fact likely to be a major mechanism for rapid 
state change in objects like GRS 1915-1-105 (Chakrabarti & Manickam 2000), is untenable in 
SS 433 because the latter is a wind accretor: thus no significant Keplerian disk is expected 
in this system to supply the soft photons, and indeed none has been detected so far (Watson 
et al. 1986; Yuan et al. 1995). 
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2.2. Resonance oscillation of accretion shocks in the presence of 

bremsstrahlung cooling 

Numerical simulations of accretion flow show that in cases where the cooling time- 
scale nearly matches the infall time-scale, a shock forms, but it then starts oscillating and 
ejects matter quasi-periodically (Langer, Chanmugam & Shaviv 1983; Molteni, Sponholz 
& Chakrabarti 1996, hereafter MSG; see Fig. lb). In order to have an oscillation period 
of around 50 s, the shock must be located at the large distance of rg^usc ~ 6400rp for a 
black hole of mass M = IOMq, where = 2GM /c^. The mass of the post-shock region 
is computed by equating the bremsstrahlung (which we assume to be the major cooling 
mechanism) cooling time and the infall time in the post-shock region (MSG): 

S Vf \ Tg J C 

where, £ is the specific thermal energy, Vfis the infall velocity and R— (7-|-l)/(7 — 1) ~4 — 7 
(these limits are for a strong shock with 7 = 5/3 and 7 = 4/3, respectively) is the compression 
ratio at the shock. Assuming the gas density (n) and temperature (T) scale as n ~ r~^/^ and 
T ~ respectively, the mass of the sub-Keplerian region of r < turns out to be 7 x lO^^g 
(with M = IOM0, 7 = 5/3). This is indeed of the same order as the mass of the bullets 
observed in SS 433. However, one has to have both the angular momentum and energy of 
the injected material comparable to the marginally bound values determined by the central 
object in order to achieve such an oscillation. On the other hand, if the mass expulsion 
from the system takes place at the similar radius of Vex ~ lO^r^ = 2 x 10^ M/Mq cm 
due to the centrifugal force, the specific angular momentum of the flow is approximately 
[rex/(2rp)]^/^rpC ~ 70rgC, which is very large compared to the marginally bound value of 
2rgC. So it is unlikely that this mechanism works in SS 433. 

2.3. Non-steady and non-linear shock oscillation 

A standing shock can form in a sub-Keplerian flow only if there are two saddle type sonic 
points and the Rankine-Hugoniot relation is satisfied at least at one point in between these 
two sonic points. However, Ghakrabarti (1990) showed that there is a large region of the 
parameter space where there are two saddle type sonic points but the shock-conditions are 
not satisfied. Even an initially supersonic accretion (such as the wind from the companion) 
can fall into this category. 

What will happen to such a realistic fiow, especially when the specific entropy at the 
inner sonic point is greater than that at the outer sonic point? RGM discovered that a fiow 
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injected with these parameters exhibits yet another type of shock-oscillation (Fig. Ic). Here 
the shock searches for a stable location and oscillates without finding it. In the first half 
of the cycle, the shock recedes far away, the post-shock region fills up, but the accretion is 
essentially completely blocked. In the second half of the cycle, the shock pushes the matter 
into the black hole, thereby evacuating the post-shock region. In a realistic simulation, 
RCM find that while the ratio of actually accreted matter to the amount available from 
the companion, = Macc/-^mj would be around 0.2 during the first half-cycle, i?ai ~ 1.3 
in the second half-cycle. The outfiow was also found to be very large. The time scale of 
oscillation was found to be Trcm ~ (4000 — 6000)rg/c for a ~ 20rg whose infall time 
is only about Tmsc — {R'^s/fgY^'^i'i'g/c) ~ (350 — 400)rp/c. Thus, this type of oscillation 
takes about a factor of Rt = T-rcm/Tmsc ~ 15 times longer than the resonance oscillation 
discussed in §2.2. For a 50 s oscillation, the location of the shock should be obtained from 
{rs,RCM/rgf^^ - (l/^)(50 s/i?r)(c/r,) ~ 10^ which gives, r,,RCM ~ 4:50rg for a IOMq black 
hole, a more physically reasonable value. Even though the size of the oscillating region goes 
down by a factor of 10 or so, compared to that involved in the resonance oscillation, the 
ejected mass need not go down (even for the same accretion rate as in the earlier case). This 
is because nearly all of the accretion fiow is accumulated in half the cycle (~ 25 s in this 
case) before being ejected (see Fig. 2 of RCM). 

Another advantage of this type of non-steady shock oscillation is that it is driven by 
centrifugal force and not by thermal cooling. Hence the result is generally independent of 
the accretion rate. Thus as long as the viscosity remains low, equivalent to having the 
Shakura-Sunyaev (1973) parameter, a < etc ~ 0.015 (Chakrabarti 1990), and Minj remains 
fairly constant, this oscillation, once established, could be sustained indefinitely. 



2.4. Magnetic rubber-band effect 

In the event of increase in magnetic activity of the disk, as could happen for instance 
when the accretion disk bends towards the binary companion during its precessional motion, 
it is not unlikely that a strong magnetic field will be first intercepted, and then advected, 
toward the inner edge of the disk. In this case the field will preferentially become toroidal 
due to shear in the rotating flow. Then, as has already been pointed out (Chakrabarti & 
D'Silva 1994; Nandi et al. 2001) the acceleration due to magnetic tension 

gt = ~ (2) 

47rr(pe + Pi) ^nrpe ' 

would be the dominant force in the post-shock region of the sub-Keplerian fiow (Fig. Id). 
Here, r is the major radius of the toroidal fiux tube and pi and Pe are the densities of the 
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medium internal and external to the flux tube respectively. The last step in Eqn. (2) is 
written because pi « Pe for a strong flux tube. Since oc 1/r and Pe oc r~^/^, we get 

ot oc r-^/^ (3) 

thus increasing rapidly as the tube comes closer to the black hole, and even surpassing the 
magnetic buoyancy, 



'Kep 



1+XL r3 J r3 ' ^'^^ 

where X = Pi/ Pe — ^ 0, and Axcp and A are respectively the specific angular momentum of a 
Keplerian disk and the disk under consideration, the accelerations in Eqns. (3) and (4) do 
cross over, since when very close to a black hole, A Axep for a sub-Keplerian flow. 

The effect of magnetic tension is dramatic, and the inner part of the disk is evacuated 
in the Alfven time scale: t/va ~ {T/aTY^"^ ^ 0.1s, for a lOM© black hole with a realistic 
Alfven speed, va — 0.1c (Nandi et al. 2001). The enhanced plasma ejection along the 
axis presumably causes sporadic magnetic flare events which would be observable as radio 
outbursts, at the same time reducing the X-ray emission from the disk which forms the 
base of the jet. Recently, such effects may have been seen (Safl-Harb & Kotani 2002) where 
simultaneous observations of 2 GHz radio and 2-20 keV X-ray ffuxes from SS 433 have been 
made, and a clear dip in X-ray ffux is seen at the same time a strong radio ffare is observed. 
It is worth noting that similar anti-correlated variations are common during flares in GRS 
1915 + 105 (Feroci et al. 1999; Naik et al. 2001) and we suggest that the ffares in SS 433 
originate in the same way. 



3. Concluding remcirks 

In this Letter, we have studied various competing processes for the creation of bullets 
which move ballistically in the jet of SS 433. We showed that blobs may be separated by: 
(1) Comptonization; (2) shock oscillations due to resonance; (3) oscillations due to inherent 
unsteady accretion solutions; (4) intense magnetic tension of the toroidal ffux tubes. We 
reject the flrst possibility because it requires a large Keplerian disk, which is unlikely. We 
are unable to distinguish at this stage which type of shock oscillation is more capable of 
producing bullet formation in SS 433, but we prefer the third possibility due to its impulsive 
and generic nature and smaller involved region. We believe that the fourth possibility of the 
inner disk evacuation should produce flaring events, but will occur rather rarely, perhaps only 
once in a single precession period, when the magnetic field of the companion is preferentially 
tilted towards the accretion disk during precessional motion. This fourth mechanism gives 
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rise to an anti-correlation between radio and X-rays, perhaps already observed in SS 433 
(Safi-Harb & Kotani 2002). 
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Fig. 1. — Four scenarios of bullet separation in SS433 are schematically shown, (a) The 
base of the jet is cooled down by soft photons from a Keplerian disk and detaches when it 
becomes supersonic, (b) Resonance oscillation of the sub-Keplerian region due to the near 
matching of the infall time with the cooling time produces discrete ejecta during the phase 
when the centrifugal barrier contracts, (c) Non-steady motion of the centrifugal barrier due 
to the inability of the flow to find a steady shock solution, (d) Magnetic tension from toroidal 
fiux tubes (shown as shaded narrow tori) causes them to collapse catastrophically in a hot 
ambient medium in rapid succession which evacuates the centrifugal barrier. The recurrence 
time of: (a) is the viscous time scale in the inner part of the disk, ~ 10s; (b-c) is ~ 50s; (d) 
is random and dictated by the enhanced magnetic activity. 



